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■ We explore progress in understanding the behaviour of cation conducting glasses, within the 
Q . context of an evolving "dynamic structure model" (DSM). This behaviour includes: in single 
, O, ! cation glasses a strong dependence of ion mobility on concentration, and in mixed cation glasses 

] a range of anomalies known collectively as the mixed alkali effect. We argue that this rich 
^ ' phenomenology arises from the emergence during cooling of a well-defined structure in glass 

• melts resulting from the interplay of chemical interactions and thermally driven ionic motions. 
^ . The new DSM proposes the existence of a new site relaxation process, involving the shrinkage 

' of empty A sites (thus tailored to the needs of A"*" ions), and the concurrent emergence of empty 
O C" sites, which interrupt the conduction pathways. This reduction of A sites is responsible in the 
^ molten glass for the sharp fall in conductivity as temperature drops towards Tg. The C sites 



X 



play an important role also in the mixed alkali effect, especially in regard to the pronounced 
asymmetries in diffusion behaviour of dissimilar cations. 



5 ' 1 Introduction 
o 
o 

Since 1990, there has been growing interest in the mechanisms of ion transport in disordered 
materials including glasses. This interest has been sustained not only by a rich diversity of 
behaviour PI I2l E] ? but also by rapid advances in spectroscopy, notably in magic angle spinning 
5^ I NMR and in EXAFS [HE], and in the increasing sophistication of molecular dynamics (MD) 
and other computational techniques 13 El El • 

Much effort has gone into accounting for "universal" aspects of glass behaviour which include 
the scaling properties of conductivity spectra [101 (HI. Thus, through the work of Dyre and 
of Funke et al |31 1121, conductivity spectra for disordered materials are well described, even 
though there is still active discussion of the underlying mechanisms. 

The important challenge, which needs to be addressed, lies in creating a conceptual framework 
to enable the behaviour of different glasses to be compared and if possible, predicted from first 
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principles. Some years ago, a start was made in this direction by Bunde, Ingram and Maass 
through introduction of the dynamic structure model (DSM). The essential idea |13|IT^n31 ITB]. 
which was based originally on EXAFS data of Greaves et al [5], is that the mobile ions play 
a decisive role in creating the sites which they occupy (as well as the doorways to these sites) 
within the glassy matrix. In essence, the mobile ions shape their immediate environments to 
meet their own requirements and also leave behind empty sites (called A for cations, B 
for i^^-cations) which act as "stepping stones" and thus define the conduction pathways. Any 
empty sites away from the conduction pathways, or where any memories of previous cations 
have been forgotten, are referred to as C sites. 

The growth of such conduction pathways explained the steep rise in conductivity in single 
ionic glasses which accompanies increases in mobile ion concentration, c, in many systems as 
a percolation phenomenon. Thus, an effective power-law behaviour often applies, where the 
conductivity scales as 

a ~ c"/'^^. (1) 

T is the temperature, k is the Boltzmann constant and a is an empirically determined parameter 
with the units of energy. This expression was deduced first from Monte Carlo (MC) simulations 
and then confirmed from an extensive survey of published experimental data ^H]- A similar 
steep increase in conductivity with ionic content also emerged from the calculations based on 
the counter-ion model |T7] and on energy landscape models ([HI, see also [Hj). 
In mixed cation glasses, discussion in terms of the DSM drew attention to the consequences 
of the strong competition between dissimilar mobile ions in their attempts to establish their 
conduction pathways. MC simulations showed how this competition led to the fragmentation of 
pathways, and hence to the sharp falls in conductivity that are among familiar characteristics 
of the mixed alkali effect (MAE). The underlying cause was attributed to dissimilar ions {A~^ 
and B^) being unwilling to visit each other's [B and A) sites, because of a mismatch between 
the requirements of the ion with what the "wrong" target site and the doorway to it could offer 
in terms of cavity size and perhaps numbers of nearest neighbours. 

The DSM continues to influence research in glassy ionics. Thus, although there are some re- 
searchers, see for example [201; who discuss the origin of mixed cation effects without invoking 
site recognition or percolation effects, there is a growing consensus that this is the right ap- 
proach. Thus, evidence for the existence of distinguishable A and B sites comes from both MD 
and molecular orbital simulations 1221 12n] • Recently, Swenson et al [21j have shown by 
reverse Monte Carlo (RMC) and bond-valence calculations how Li~^ and Rb'^ ions in mixed- 
cation phosphate glasses find themselves on pathways containing mixtures of Li and Rb sites. 
They too attribute the loss of ion mobility to a mismatch energy along the lines previously 
proposed by some of us [T3] . 

As a consequence of this activity, new questions continue to arise. First and foremost, there is 
the role of the network. How much does the matrix participate in inter-site conversions {A to 
C, etc.)? Do these processes occur above and below Tg, or only above Tg as suggested recently 
by Maass How closely linked are the chemical and physical processes which lie behind the 
behaviours observed in molten and glassy materials? Are we yet in a position to say that we 
have a theory of the mixed alkali effect? The present paper is an attempt to see how far these 
questions can be answered in the light of an updated dynamic structure model. 
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Figure 1: Ionic conductivity as a function of temperature for (a) Na20 — 2Si02 and K2O — 2Si02 and 
(b) 0.7 Agl — 0.3Ag2MoO4 glasses in both the molten and the glassy state. (Redrawn from Ref. j26j). 



2 Single cation systems 

2.1 Different temperature regimes 

Figure n contains Arrhenius plots showing how the conductivity varies exponentially with re- 
ciprocal temperature for Na20 — 2Si02, K2O — 2Si02 and 0.7AgI — 0.3Ag2MoO4 glasses in 
both the molten and glassy states j^EI- Characteristically, the plots are curved above Tg and 
linear below. For convenience, we identify three regions of interest. First, there is the high tem- 
perature regime, where conductivities approach those of typical molten salts (ca. lQ~^cm~^). 
Second, there is a region of under cooling which extends right down to Tg. This region is asso- 
ciated with huge increases in melt viscosity (from about 10^ to 10^'^ Pa.s) and by comparison 
a rather smaller decrease (only two to three orders of magnitude), in the ionic conductivity. 
Third, there is the region below Tg in which ion transport proceeds with a constant activation 
within an apparently rigid solid. The conductivity decreases comparatively rapidly with 1/T 
in the middle region, while it decreases comparably slower in the first and third regimes. 
At temperatures far above Tg, (region 1), all ionic motions are coupled to processes involved in 
viscous flow. These microscopic processes include the switching of non-bridging oxygens (nBOs) 
between neighbouring silicon atoms and temporary changes in the coordination number of Si 
from 4 to 5 or even to 6 ^7\. The sites created by the mobile ions may be of limited duration, 
but they will be precisely determined by the interplay of forces representing the chemical 
requirements of both mobile ions and the various silicate or other anionic species present in the 
melt. 

However, it is known from the work of Moynihan and Angell ^Hl 1^ that undercooling leads 
to progressive decoupling of cation motions from the host matrix, with an associated increase 
in the ratio of relaxation times for viscous flow to those for conductivity, i.e. of Tghoar/Tsigma- 
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The value of this ratio at Tg is called the decoupling index, R^., and reaches values around 
10^^ in typical silicate glasses. In practice, this means that since structural relaxation times 
are approximately 100 s at Tg, (as measured, for example, in a typical differential scanning 
calorimetry (DSC) experiment), the corresponding electrical relaxation time (r = RC) is about 
10"^ s. In simple terms, this gives an indication of how long ions are remaining in their sites 
at the glass transition temperature. The inverse of r is the average effective hopping rate. 
On further cooling below Tg, the structure changes much less, but ions spend increasingly longer 
times in their sites. At room temperature, a sodium ion in a typical silicate glass might sit for 
approximately 1 s in its site before moving on. The corresponding relaxation time for viscous 
flow is so long it would better be expressed on geological time scales [20] • 
In the present context, it is noteworthy that Binder et al also draw attention to the 
second temperature region described above, which they identify as falling between the critical 
temperature Tc of mode coupling theory (MCT) where the system undergoes some "structural 
arrest" and the "solidification temperature" (i.e. the experimental Tg). They also point out 
that the time scale for viscous flow at is about ten orders of magnitude shorter than it is at 
Tg, and suggest that it is unlikely that significant changes will not occur as a liquid is cooled 
between and Tg. This is precisely the region of melt behaviour where the decoupling ratio 
is increasing from values close to unity to much larger values. Indeed, within the framework of 
MCT it seems likely that near T^ the dynamics of the network forming species will be dominated 
by mutual blocking effects, and their diffusion will be drastically reduced. In effect, the network 
structure is forming around T^, and already a pattern of behaviour is being established which 
continues on cooling even into the solid glass. 

2.2 The mechanism of ion transport 

Evidence from recent NMR studies jH2] shows that cations in glass will typically find sites near 
at least one nBO. More generally, in glasses where nBOs may be absent, e.g. in aluminosilicate 
glasses, the cations will be found in regions close to where negative charge resides in the network. 
The increasing values of Tshear/Tsigma, see above, show that even in the melt the local network 
structure persists for much longer times than the ions reside in their sites. We are looking at 
a hopping process which is well established in the molten state. 

But what kinds of site do the ions hop into? What distinguishes an A from a C site has not 
so far been rigorously defined. It was assumed [15J that A character expresses a "goodness of 
fit" , and combines a number of factors that determine the accessibility of the site to an A"*" ion, 
including cavity size, coordination number, and the availability of negative charge as indicated 
often by the presence of nBOs. 

We now propose to distinguish two kinds of C sites. First, there are C sites, which in silicate 
melts would be remote from nBOs. These can only become A sites with the active assistance of 
chemical rearrangements occurring in the silicate matrix. It is likely that such rearrangements 
occur readily only at high temperatures (possibly in the vicinity of the Tc of MCT), and that 
they occur less and less frequently as the melt is cooled towards Tg. 

The second kind of C site has very different origins. Let us suppose that a freshly vacated A 
site looks just like a filled site, but with the cation missing. As long as it stays like that, we 
shall call it an A site. However, after being vacated, the A site may start to relax in such a 
way as to minimise the (now largely uncompensated) repulsions between negatively charged 
oxygens in what was previously the primary coordination sphere of the cation. This relaxation 
might be envisaged, for the purpose of illustration, as the reorientation oi Si — O bonds to 
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Figure 2: Site relaxation, showing the collapse of an empty A site, either during cooling of molten 
glass or when glass is subjected to external pressure, with the consequent creation of a new C site 
and the resulting disappearance of some free volume at the same time. 

avoid pointing the oxygens directly at each other (as will be the case in the newly emptied A 
site) but instead towards neighbouring silicon atoms which (because of differences in chemical 
electronegativity) will carry small positive charges. Whatever the actual (microscopic) details 
of this relaxation process, we call this new site a C site. 

The processes of site relaxation occur in the molten state, but because melts are relatively rich 
in "free volume" , the C and A sites may not differ too much from each other. However, the melt 
contracts on cooling and so some free volume must disappear. We suggest that this contraction 
occurs most readily by the conversion of empty A into smaller C sites, as illustrated in Fig. |21 
Below Tg, the total number of A and C sites will remain constant and further shrinkage of the 
C sites will be less important. 

Below Tg, two alternative scenarios may in fact be envisaged. In the first, the structure may 
be completely frozen, with all sites retaining their A, C, and C identities as they were defined 
at Tg. Such a situation would be consistent with many contemporary views of glass behaviour 

[nii2si. 

The alternative would be that minor structural fluctuations such as those involved in converting 
C into A sites or vice versa, would still occur, since changes in network topology would not be 
involved. The main difference between C and A sites would thus be related to local density 
fluctuations. In this way, it would still be possible to refer to a dynamic energy landscape, 
which is attributed to the network structure, even below Tg. 

Whichever scenario is adopted, this new description has several advantages. As one example, 
it allows one to visualise what is happening in the Arrhenius plots (see Fig. ^ in different 
temperature regimes. Above Tg, the progressive increase in slope seen during cooling, is a 
consequence of the decreasing number of empty A sites as they are converted into smaller C 
sites. The arrest at Tg signals the effective end of this process. Below Tg, the energy landscape 
is determined by the distribution of empty A and C sites. 

We can go on further to consider that the number of C sites is significantly greater than 
the number of empty A sites. The situation in glass has some resemblance to the "vacancy- 
like" models proposed recently by several authors PSHSHIISII- Under the influence of external 
pressure it may even be possible to convert some of the remaining empty A into smaller C sites. 
Elsewhere, one of us [3^ argues that the apparent "activation volume" for charge transport 
in cation-conducting glasses can be interpreted in this way, and indeed is equal to the volume 
difference between empty C and A sites. 

Regarding the existence of prefered pathways or conducting channels jSHl EZl ISH] : In a recent 
paper, Meyer et al. identify a prepeak observed in both neutron scattering data and 
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Figure 3: (a) Diffusion coefficients in Na/Cs trisilicate glasses at 396 °C as a function of tlie com- 
position dependence Cs/{Cs + Na). (Data after jlTj) redrawn from Ref. PP). (b) Numerical results 
from the DSM for the normalized diffusion coefficients Da{x)/Da{^) and Db{x)/ Db{0) of cations A 
and B as a function of x for a symmetric choice of parameters and for three different temperatures. 
The full symbols refer to the A ions, the empty symbols to the B ions (redrawn after Ref. jl5j). 

molecular dynamics simulations with the emergence of Na"^ ion channels in sodium silicate 
melts and glasses. In terms of the new DSM, we expect that Na and C sites close to the nBOs 
form these channels. Hence, suppression of the formation of nBOs, e.g. by addition of AI2O3 
to sodium silicate glasses should disrupt the pathway structure, in agreement with very recent 
experiments by Kargl et al. [SHI- 

3 Mixed cation glasses 

3.1 The mixed alkali effect (MAE) 

The mixed alkali effect HHl 1^ embraces a wide range of phenomena. Some of these, 
such as diffusivity crossovers and conductivity minima are directly related to changes in ion 
mobility, and are associated with maxima in activation energies and entropies, and in the 
activation volume |12I. There are, however, other anomalies, which only indirectly involve 
ion transport. These include minima in Tg jlHl HHj; large internal friction peaks [101 IS] and 
volumetric relaxations [101 133] observed below Tg, and viscosity minima observed in molten 
glass lini- 

Figure IHl^a) shows the diffusivity crossover in Na — Cs silicate glasses reported by Jain and 
coworkers • The behaviour of Na differs markedly from that of Cs. While the self-diffusion 
coefficient of the latter decreases exponentially with decreasing Cs content, the corresponding 
coefficient for Na levels off at low Na contents. This asymmetry is disguised in conductivity 
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plots, which tend to be rather U or V-shaped in appearance. The conductivity is, however, 
always dominated by the more mobile and numerous (majority) cations; the significant differ- 
ences in behaviour between large and small ions tend to show up more clearly only when the 
ions in question are in the minority. 

On the other hand it is possible to overstate the degree of asymmetry. Measurements of 
activation volumes {Va = —RTdlna/dP, with the pressure P and the gas constant R) in 
mixed cation glasses reveal increases in Va occurring at both ends of the composition range, 
when larger ions are diluted by smaller ones or vice versa Ell- This evidence for some 
degree of symmetry has suggested some degree of coupling between the local motions of the 
ions, involving ions moving into B sites and B'^ ions moving into A sites j42| l4Hj . 
In the dynamic structure model, the asymmetry was expressed in terms of different jump rates 
of A'^ and B~^ cations. In the published simulations for simplicity and for revealing the 
crucial effect of mismatch between unlike ions, a hypothetical case was considered of a glass 
containing a mixture of cations, A'^ and B^, alike in all respects except that a mismatch 
energy appeared whenever either ion entered sites belonging to the other. Figure Efb) shows 
the crossover in diffusivities predicted by the quantitative model ^H]- This figure correctly 
reveals the steep falls in ion mobility found in the dilute foreign alkali region, but of course it 
could not replicate the contrasting behaviours of Na and Cs shown in Figure a). 
It is now appropriate to look again at the MAE in the light of the newly " updated" dynamic 
structure model. 

3.2 The MAE in molten glass 

There have been several reports in the literature of the MAE found in molten glass [101 HO] • 
Baucke and Werner have shown, that for a series of Na20 — K2O — CaO — Si02 glasses, there 
are significant deviations from additivity in ionic conductivity, and especially in the activation 
energy well above Tg (see Fig. |H^a) and (b), respectively). The effect is thus still apparent 
in melts of relatively modest viscosity, arguably close to the critical temperature, Tc, of mode 
coupling theory. The persistence of the MAE up to such high temperatures is strong evidence 
that distinctions between A and B sites can still be made. 

According to our updated dynamic structure model, the sites are formed in the melt between 
Tc and Tg where the volumes of the A and the C sites do not differ very much. In mixed 
alkali melts, the situation is further complicated owing to the different sizes of the A'^ and 
the B^ ions. When the smaller A'^ ions are in the overwhelming majority, the C sites are on 
average smaller than when the larger B~^ ions are in the majority. Close to the conduction 
minimum, where nearly the same numbers of and B^ ions are present, we would expect 
a broad distribution of C sites, including sites to be found in each of the single alkali melts. 
Under such conditions, ion transport involves continual melt rearrangement, and this gives rise 
to the maximum in activation energy which is so evident in Baucke and Werner's data. 

3.3 The MAE in solid glasses 

The starting point, as in earlier discussions |15j, is to recognise the mismatch effect which arises 
whenever ions enter "wrong sites". It is clear that larger {B~^) ions will find it harder (as a 
natural consequence of electron repulsions) to enter smaller (A) sites than vice versa. Also, 
because of the broad distribution of C sites which formed in the melt, there are many more 
C sites available to A^ ions than there are to B^ ions. Accordingly, in Na/Cs glasses (see 
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Figure 4: (a) Conductivity as a function of the potassium mole fraction x in the molten mixed- 
alkali system (1 — x)Na20 ■ XK2O ■ 0.7CaO ■ 4:.8Si02 at various temperatures, (b) The corresponding 
activation energies at various temperatures taken from curved Arrhenius plots. (Redrawn after iBOj ) . 

Fig-E^a)), Cs ions will be expected to suffer much larger losses in mobility than Na ions, which 
indeed is the experimental result. This explains also why Na ion mobilities level off at low Na 
contents, but to a value lower than that of the majority Cs ions (if ions jump into sites which 
are too large for them, they will tend to "roll back" into their original sites.) 
If cations do not differ too much in size, as in Li/Na and Na/K mixtures for example, the A'^ 
cations may find themselves entering B sites in preference to less attractive C sites. In this 
way, ions and A sites become randomly mixed up with ions and B sites (as indeed is 
indicated both by infrared and NMR spectroscopies) jH |^ . The frequent entry of A'^ ions into 
B sites and vice versa leads to a loosening of the structure. This loosening was seen by Ingram 
and Roling as being the underlying cause of minima in Tg and of melt viscosity jlH]; as 
interpreted with the concept of matrix mediated coupling (CMMC). 

Remarkably, these indirect effects may be smaller when cations differ more markedly from 
each other. This was found in internal friction experiments [Ulj. Very small loss peaks were 
reported in Li/Cs silicate glasses, compared with much larger ones in Li/Na glasses. We now 
understand this additional anomaly (i.e. the suppression of an important aspect of the mixed 
alkali effect) in terms of Li~^ and Cs+ ion simply differing too much from each other to enter 
each other's sites. 

4 Conclusions 

We introduce an updated dynamic structure model for explaining the complex behaviour of 
single and mixed cation glasses. As in the original version jElEj, we distinguish between A, 
B and C sites. However, in contrast to the old DSM, we assume that A and B sites can only 
be generated close to the countercharges (which often are non-bridging oxygens). When an A 
site is vacated, it relaxes to a C site whose size is determined largely by the host matrix, but 
it is better suited to accommodate mobile ions than are the C sites, which are remote from 
the localized countercharges. In a single cation glass, the A and C sites found close to the 
countercharges, form "stepping stones" for the mobile A'^ ions. 

As in the recent papers j^Hl EH EH]; we expect that these pathways will appear in the melt, 
perhaps at a temperature close to of mode coupling theory, when the diffusion of the network 



8 



forming units is blocked, and viscosity becomes very temperature dependent. These pathways 
constitute an integral part of glass structure below Tg. The pathways become interrupted when 
the stepping stones are too far apart, i.e. when the contercharge concentrations are too small. 
The pathways may be "frozen in" at Tg or they may be allowed to fluctuate in the glassy state 
as a result of localised density fluctuations. In both cases, the loss of good pathways leads to 
the abnormal decreases in conductivity with decreasing ion content, an effect which is especially 
severe in certain systems such as alkali borate glasses. 

In a mixed ionic glass, it is much more difficult to find separate pathways for the A"^ and the 
ions, made up of A and C and of B and C sites, respectively. This appears to be the underlying 
cause of many of the observed anomalies. As in the old DSM, we envisage that there will be a 
mismatch whenever a (small) A~^ ion enters a (large) B site or vice versa. We now recognise 
that this mismatch is greater if a larger cation is trying to enter a smaller site, and in this 
way we can account for the observed asymmetries in diffusion behaviour, these being especially 
significant with cations of markedly different sizes, as in the Na — Cs system for example. By 
contrast, the smaller of the two ions {A'^) will make use of available C sites, when the supply 
of {A) sites is insufficient to provide extended diffusion pathways. This comment is consistent 
with the "levelling off" in the diffusion coefficient of Na^ ions at low concentrations, as seen 
in the Na — Cs system. 

However, in molten glass, the ^4+ and B^ ions will have sufficient energy to visit each other's 
sites, and regularly disturb the local structure. This could possibly explain the strong reduction 
of Tg and of melt viscosity reported in many mixed alkali systems. When the ions do not differ 
too much in size, these site exchange effects may persist below Tg. In this way, we could account 
for the large internal friction peaks seen in Li/Na silicate glasses and for the anomalous increases 
in activation volume which point to some coupling between the motions of the A^ and the B^ 
cations. 

What is urgently needed is more information about empty sites in glass, and how the nature 
of these sites depends on glass composition and temperature. In the related field of polymer 
electrolytes, diagnostic information is coming from positron annihilation lifetime spectroscopy 
or PALS ^52,,53j. This technique exploits the variable lifetimes of the ortho-positronium species 
in order to determine the size of available cavities within the host matrix. Such information, to- 
gether with activation volumes determined from variable pressure measurements, might enable 
A and B sites to be distinguished from each other and from the possibly more numerous C 
sites. Comparisons made between fast quenched and slow cooled glasses would provide infor- 
mation not only on structural dynamics but also on the growth and consolidation of percolation 
pathways. 

Finally, it is interesting to note previous comments in the literature to the effect that the 
concentration dependence of conductivity in mixed cation glasses is closely related to those 
encountered in single cation systems [Hlj. If this rule applies quite generally, it would imply 
that the C sites encountered in single cation glasses must be very similar to those encountered 
in mixed cation systems. This may or may not be the case. Further consideration of the MAE in 
such glass forming systems as the aluminosilicates [531 EH] , where the conductivity dependence 
in single cation glasses is weak [HZj, should be highly informative. 
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